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Abstract 

The unique existence of a weak solution to the homogeneous closed 
Dirichlet problem on certain D-star-shaped domains is proven for a 
mixed elliptic-hyperbolic equation. Equations of this kind arise in 
models for electromagnetic wave propagation in cold plasma. A re- 
lated class of open boundary value problems is shown to possess strong 
solutions. 
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1 Introduction 

Boundary value problems for mixed elliptic-hyperbolic equations may be ei- 
ther open or closed. In the former case, data are prescribed on a proper 
subset of the boundary whereas in the latter case, data are prescribed on 
the entire boundary. It is shown in Sec. 3 of [22J that if k = 1/2, the closed 
Dirichlet problem is over-determined for the equation 



x 



UXX + Uyy + KU X = (LI) 



on a typical domain, where u (x, y) is required to be twice-continuously differ- 
entiable on the domain. However, this equation arises in a qualitative model 
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for electromagnetic wave propagation in an idealized cold plasma ([32], eq. 81; 
see also [30J, eq. (9)). Physical reasoning suggests that the closed Dirichlet 
problem for (II. ip should be well-posed in a suitable function space, at least 
for some choice of lower order terms; so the result reported in [22] would 
appear to represent a serious deficiency in the physical model. See Sec. 1 of 
[22] for a discussion, in which the problem of formulating a closed Dirichlet- 
like problem that is well-posed in an appropriate sense is characterized as an 
"outstanding and significant problem for the cold plasma model." 

Except for the point at the origin, eq. ( II. ip can be transformed into an 
equation of Tricomi type; see, e.g., Sec. 2.4.3 of |26j. Such equations have 
somewhat more desirable analytic properties than eq. (II. ip . However, both 
the physical and mathematical interest of eq. (II. ip arise from the tangency 
of the resonance curve to a flux line at the origin. This is the point at which 
plasma heating might occur in the physical model, and a point which appears 
to be singular in studies of the solutions to (jl.ip ; see [22], [30J, and [39]. Thus 
it is not generally useful to transform eq. (Il.ip into an equation of Tricomi 
type. 

Using methods introduced by Lupo, Morawetz, and Payne [T7], [T8] for 
equations of Tricomi type, we show in Sec. 2 the weak existence of a unique 
solution to a homogeneous closed Dirichlet problem for the formally self- 
adjoint (« = 1) case of eq. (II. ip . This extends a recent result [21] in which 
the existence of solutions having various degrees of smoothness was shown in 
certain cases to which uniqueness proofs did not seem to apply. At the same 
time, it extends the unique-existence arguments in [17] to an equation which 
is not of Tricomi type. 

Another well known problem in elliptic-hyperbolic theory is the determi- 
nation of natural conditions for boundary geometry; see the discussions in [TJ, 
[2D] . [21], [21], EB]) an d [2H]- Heuristic approaches to determining boundary 
geometry tend to focus on physical [19] or geometric [2S] analogies for the 
specific equation under study. In his theory of symmetric positive systems [8], 
Friedrichs proposed intrinsic mathematical criteria for the well-posedness, or 
admissibility, of boundary conditions. But Friedrichs' conditions are also tied 
to the specifics of the particular symmetric positive equation under study and 
are algebraic rather than explicitly geometric. We will require boundary arcs 
to be starlike with respect to an appropriate vector field. This approach to 
boundary geometry was introduced by Lupo and Payne [16J. We note that 
algebraic conditions in certain very old results can be reinterpreted as the 
requirement of a starlike boundary; see, e.g., [UJ. Our results provide further 
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evidence that domain boundaries which are starlike in this generalized sense 
are natural for elliptic-hyperbolic boundary value problems. 

In Sec. 3 we investigate the solvability of open boundary value prob- 
lems for a class of symmetric positive systems on domains having starlike 
boundaries. This class includes an equation originally studied by Cinquini- 
Cibrario in the early 1930s. The equation arises in the so-cold "slab" model 
of zero-temperature plasma and in models of low-temperature atmospheric 
and space plasmas (Sec. 3.1). 

The boundary conditions in Sec. 3 are mixed in the sense that a Dirich- 
let condition is placed on part of the boundary and a Neumann condition 
is placed on another part. However, our methods also apply to the case in 
which either a Dirichlet or a Neumann condition is imposed over the entire 
elliptic boundary. Because the boundary value problems considered in Sec. 
3 are open, the results of that section may be less interesting physically than 
those of Sec. 2. But open boundary conditions can be expected to imply more 
smoothness on the part of solutions than is obtained from closed boundary 
conditions, and we provide conditions for the existence of solutions which 
are strong in the sense of Friedrichs. Conditions for the existence of strong 
solutions to elliptic-hyperbolic boundary value problems were also discussed 
in Sec. 3 of [24], but briefly and inadequately. Section 3 of this report revises 
and extends (to the open case of mixed and Neumann problems) the treat- 
ment of strong solutions in [24] . The existence question for weak solutions to 
open boundary value problems for equations of the form (11. ip was considered 
in [23] and @T]. 

1.1 Remarks on the physical model 

In the cold plasma model, the plasma temperature is assumed to be zero 
in order to neglect the fluid properties of the medium, which is treated as 
a linear dielectric. In the case of wave propagation through an underlying 
static medium having axisymmetric geometry, equations of the form (11. ip 
model the tangency of a flux surface to a resonance surface. At the point of 
tangency, plasma heating might occur even in the cold plasma model [3H] • I n 
two dimensions, flux surfaces (level sets of the magnetic flux function) can 
be represented by the lines x = const., and a resonance surface (frequencies 
at which the field equations change from elliptic to hyperbolic type) by the 
curve x = y 2 . In such cases a plasma heating zone could lie at the origin of 
coordinates. This conjecture is supported by numerical [22] and classical [30] 



3 



analysis which suggests that the origin is a singular point of eq. (11. ip . 

The main physical implication of our result is that, although the Dirichlet 
problem for the cold plasma model is manifestly ill-posed in the classical 
sense, there is a weak sense in which the closed Dirichlet problem is no longer 
over-determined: a unique weak solution to the closed Dirichlet problem 
exists in an appropriately weighted function space. The weight function 
vanishes on the resonance curve. In Corollary 6 of [21] an existence theorem, 
without uniqueness, was demonstrated in a weighted function space in which 
the weight function vanished on the line y = 0. Thus in each case the weight 
function vanishes at the origin of coordinates. These results support the 
physical conjecture of a heating zone at that point. 

For discussions of the physical context of eq. ( 11. ip . see Sec. V of [39J and 
Sec. 4 of [10] , in which eq. ( II -ip with k = is proposed as a qualitative model 
for erratic heating effects by lower hybrid waves in the plasma. See also [30J, 
in which a model for electrostatic waves in a cold anisotropic plasma with a 
two-dimensional inhomogeneity yields, by a formal derivation, an equation 
for the field potential which is similar to (II. ip . Briefly, the derivation pro- 
ceeds as follows: at zero temperature, the field equations reduce to Maxwell's 
equations, which are written for the electric displacement vector D. The com- 
ponents of D are written in terms of the dielectric tensor for the medium. 
The coordinate system is chosen so that, in two dimensions, the y-axis is 
collinear to a longitudinally applied magnetic field. In that case the dielec- 
tric tensor assumes a particular form which leads, when the resonance curve 
is tangent to the flux line at the origin, to an equation having the form con- 
sidered here. Precisely, the equation derived in [30] is eq. f !3.29|) of Remark 
in), Sec. 3, below, with particular choices of «i, «2, &(y), Ki, and k 2 . 

The derivation of the governing equation for the case of fully electromag- 
netic waves is much more lengthy, and occupies much of the memoir [3H] ■ In 
that case, equations having the form (II. ip provide only a qualitative descrip- 
tion of the physical model. Section 2.5 of [26] is devoted to a brief derivation 
of eq. (11.11) for electromagnetic waves. 

The classic text for wave propagation in plasma is [55] ; see also [9]. See 
Ch. 2 of [36] for a concise introduction to the physics of electromagnetic 
waves in cold plasma, and [26] for a recent review of mathematical aspects 
of the cold plasma model. 

In the sequel we assume that fl is an open, bounded, connected domain 
of 1R 2 having at least piecewise continuously differentiable boundary with 
counterclockwise orientation. For nontriviality we require that Q contain an 
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arc of the resonance curve x = y 2 (x = a(y) in Sec. 3). Additional conditions 
will be placed on the domain where required. 



2 Weak solutions to closed boundary value 
problems 

Following Sec. 3 of [17] we define, for a given C 1 function K (x, y) , the space 
L 2 (fl; \K\) and its dual. These spaces consist, respectively, of functions u for 
which the norm 



\ u \\L 2 (n;\K\) 




\K\u 2 dxdy 



1/2 



is finite, and functions u G L 2 (Q) for which the norm 



l^Q;!^!- 1 ) 




K\ l u 2 dxdy 



1/2 



is finite. Analogously, we define the space K) to be the completion of 

C^°(fi) with respect to the norm 



\U\\m(Q,;K) 



J J (\K\u 2 x + u 2 + u 2 ) dxdy 



1/2 



(2.1) 



and introduce the space Hq(Q; K) as a closure in this space. Using a weighted 
Poincare inequality to absorb the zeroth-order term, we write the Hq(Q; K)- 
norm in the form 



U Ml 



Hl(p.;K) 



J Jj\K\u 2 x + u 2 y ) dxdy 



(2.2) 



The dual space H 1 (Q; K) is defined via the negative norm 

\(w t (p)\ 



\W\\H-i(n;K) = SUp 

o^ecg°(n) 



where ( , ) is the Lax duality bracket. 

Various lower-order terms have been associated to eq. (11. ip in the liter- 
ature on the cold plasma model; only the higher-order terms contribute to 
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resonance. These lower-order terms were considered explicitly in Sec. 2 of 
[23], in which we showed the existence of L 2 solutions to a closed Dirichlet 
problem for eq. (jl.ip . Here we consider only the formally self-adjoint case 

Lu = [K (x, y) u x ] x + u yy = f (x, y) , (2.3) 

for the type-change function K (x,y) = x — y 2 . As in [23], the inhomogeneous 
term / (x, y) is assumed known. For this special case of lower-order terms, 
we are able to show the existence of solutions to the closed Dirichlet problem 
having much higher regularity than was the case in [23]. 

The forcing function / in eq. (12.31) (and also in eq. (13.1 1) , below) could arise 
physically under appropriate conditions on the charge density. In addition, 
the presence of / provides mathematical generality, avoids trivial solutions, 
and allows the conversion of solutions to the Dirichlet problem with homoge- 
neous boundary conditions into solutions to a problem for the homogeneous 
equation with inhomogeneous boundary conditions. Finally, analysis of the 
inhomogeneous equation is useful when, as in our case, one expects singular- 
ities in the homogeneous equation. 

In accordance with standard terminology, we will often refer to the curve 
K = on which eq. (I2.3P changes type as the sonic curve. This terminology 
is borrowed from fluid dynamics; in the context of the cold plasma model, 
the sonic transition occurs at a resonance frequency. 

Following Lupo, Morawetz, and Payne [17], we define a weak solution of 
eq. (12. 3p on Q, with boundary condition 

u(x,y) =0\/(x,y) Edtt, (2.4) 

to be a function u G Hq(Q] K) such that V£ G Hq{VL; K) we have 

(Lu, £) = - J J {Ku x £ x + u y Q dxdy = (f, f ), 

where ( , ) is duality pairing between Hq (Q, K) and H~ l (Q, K) . In this case 
the existence of a weak solution is equivalent to the existence of a sequence 
u n G C™(Q) such that 

I \ u n ~ u\\ H i {n . K ) -)> and \ \Lu n - f\\H-l(Cl;K) -> 

as n tends to infinity. 
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Following Lupo and Payne (Sec. 2 of [16J), we consider a one-parameter 
family ip\ (x, y) of inhomogeneous dilations given by 



^ (x,y) = (\- a x,\-Py), 



where a, {3, A G M + , and the associated family of operators 



^ x u = u o ip x = u x - 



Denote by D the vector field 




J|A=1 



— axd x — {3yd, 



(2.5) 



An open set Q C M 2 is said to be star-shaped with respect to the flow of D 
if V (x , y ) G Q and each t G [0, oo] we have F t (x , y ) C Q, where 



If a domain is star-shaped with respect to a vector field D, then it is possible 
to "float" from any point of the domain to the origin along the flow lines of the 
vector field. If these flow lines are straight lines through the origin (a = (3) , 
then we recover the conventional notion of a star-shaped domain. By an 
appropriate translation, the origin can be replaced by any point (x s ,y s ) in 
the plane as a source of the flow. In that case we obtain a translated function 
F t for which 



Moreover, whenever a domain is star-shaped with respect to the flow of a 
vector field satisfying (12. 5p . the domain boundary will be starlike in the sense 
that 



where n is the outward-pointing normal vector on the boundary dQ. See 
Lemma 2.2 of [TB]. In equivalent notation, given a vector field V — — (b,c) 
and a boundary arc T which is starlike with respect to V, the inequality 



F(x ,y ) = {x{t),y{t)) = (rroe-^j/oe-^) . 



lim F t (xo,y ) = (x s ,y s ) W(x Q ,y ) G Q. 



(ax,(3y) ■ n(x,y) > 



brii + cn 2 > 



(2.6) 



is satisfied on T. 
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2.1 An auxiliary problem 

We employ the dual variational method, an integral variant of the abc method, 
introduced by Didenko [7J and developed by Lupo and Payne [IS]. (Note 
that the term "dual variational" also appears in elliptic variational theory, 
in which it means something entirely different.) Denote by v a solution to 
the boundary value problem 

Hv = u in Q (2.7) 
for u G Co°(fi), with v vanishing on <9f2\{(0, 0)}, 

lim v(x,y)=0, (2.8) 

and 

Hv = av + bv x + cv y . (2.9) 

Assume that Q is star-shaped with respect to the flow of the vector field 
V — — (b, c) ; b — mx and c = \iy\ fi and m are positive constants and a is 
a negative constant; the point (x,y) = (0,0) lies on dQ. Step 1 in the proof 
of Lemma 3.3 in [T7], which treats the harder case of a non-differentiable 
coefficient, demonstrates that v satisfying (I2.7p - (I2.9I) exists and lies in the 
space C° (H) n Hi (O; K) . The proof is straightforward in our case, except 
perhaps for the justification of (12. 8p . which we outline: 

Because Q is star-shaped with respect to — (6, c) , any flow line entering 
the interior of Q from dfl will remain in f2. So the method of characteristics 
yields a smooth solution to (12. 7L (12. 9ft . with a possible singularity at the ori- 
gin. It remains only to analyze the limiting behavior of the solution near the 
origin. Because the origin lies on the boundary and u has compact support, 
we can restrict our attention to an e-neighborhood N £ of the boundary, 

N £ (dn) = {(x,y)eU\dist((x,y),dn)<e}, 

where e is so small that N e lies in Q but outside the support of u in Q. In 
this subdomain we solve the Cauchy problem for the homogeneous equation 

mxv x + fiyVy = \a\v. (2-10) 

The characteristic equation for (I2.10p . given by mxdy = fiydx, can be inte- 
grated to yield 

y = cx> M/m , (2.11) 
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where c is a constant. The method of characteristics thus yields a general 
solution in N e having the form 

v(x,y) = <p0^yW'>, (2.12) 

where (p is an arbitrary C 1 function that may be prescribed along a non- 
characteristic curve in N e . Following the evolution of the solution (j2.12p along 
the curves ( 12. lip , we obtain ( 12. 8 p as in eq. (3.13) of [17] . 

Following Sec. 2 of [7J (see also the Appendix to [IS]), the dependence of 
the solution v to the problem ( |2.7p - (j2.9p on the forcing function u is repre- 
sented by an operator X, writing v = Xu. We have the integral identities 

(Xu, Lu) = (v, Lu) = (v, LHv) . (2.13) 

A good choice of the coefficients a, m, and \x in the operator H on the right- 
hand side of this identity will allow us to derive an energy inequality, which 
will be used to prove weak existence via the Riesz Representation Theorem; 
see Ch. 2 of [2] for a general treatment of such "projection" arguments. 

2.2 Main result 

The following is a small but crucial extension of [21], Theorem 5. 

Lemma 2.1. Suppose that x is non-negative on Q and that the origin of 
coordinates lies on dQ. Let Q be star-shaped with respect to the flow of the 
vector field V = — (b, c) for b = mx and c = fiy, where m and \l are positive 
constants andm exceeds 3/i. Then there exists a positive constant C for which 
the inequality 

\\u\\L 2 (n ; \K\) < C\\Lu\\ H -t(p,;K) 
holds for every u G (fl) , where K(x,y) = x — y 2 and L is defined by 

Proof. Let v satisfy eqs. (12.7p - fj2.9n on Q for a = —M, where M is a 
positive number satisfying 

2 
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for some sufficiently small positive number S. Integrate the identities ( I2.13P 
by parts, using Prop. 12 of |24J and the compact support of u. We have 

v ■ LHv dxdy = \ f {Kv 2 x + v 2 ) (cdx — bdy) 
2 Jan 

+ avl + -fv 2 y dxdy, (2.14) 
J Jn 

where 




a = K[ ^A-a] +t h + L KyC 



(t ~ ^ ~ 6 ) x + 6y2 



and 



2 2 2 ^ ^ 



On the elliptic region Q + , K > and 

^— — /i — 5 j x > (J^ — S j x > Sx 

provided we choose 5 so small that /i/4 > 5. Then on Q + , 

u> 5 (x + y 2 ) > 5 (x - y 2 ) = 5K = 5\K\. 
On the hyperbolic region Q~ , K < and 

a = (—-^x + 6(y 2 -x)>^x + 6 {-K) > 5\K\. 

We show that the integrand of the boundary integral in (I2.14p is non- 
negative: Because v G C° (H) n Hi (Q; K) , 

\K\v 2 x + v 2 y = 

in a sufficiently small neighborhood of dfl. This implies that v 2 = v 2 = in 
a sufficiently small neighborhood of dQ\ {K = 0} . That in turn implies that 

Kv 2 x + v 2 y = 

in a sufficiently small neighborhood of dQ\ {K = 0} . On the curve K = 0, 
we have x = y 2 and cfe = 2ydy, implying that 

cdx — bdy = (2fiy 2 — mx) dy = (2/i — m) xdy. 
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But m > 3/i > 2//, i > on fi, and dy < on i^; so on the resonance curve, 

^ (crfx — bdy) > 0. 

Thus the integrand of the boundary integral in (I2.14p is bounded below by 
zero. 

We find that if 5 is sufficiently small relative to //, then 

(v,LHv)>5 [ [ (\K\v 2 x + v 2 y )dxdy. (2.15) 
J Jn 

The upper estimate is immediate, as 

(v, LHv) = (v, Lu) < \\v\\ H i {Q . K) \\Lu\\ H - 1(n . K) . (2.16) 

Combining (" 12 . 1 5[) and (I2.16p . we obtain 

Ikllfl-^n;!?) - C WLuWh-x^.k} ■ (2-17) 

The assertion of Lemma I27H now follows from (12. 7p by the continuity of H as 
a map from Hq (Q; K) into L 2 (Q; \K\) . This completes the proof of Lemma 
I2TT1 

Theorem 2.2. Let Q be star-shaped with respect to the flow of the vector field 
—V = (mx,fiy) , where m and /i are defined as in Lemma \2.1\ Suppose that 
x is nonnegative on Q and that the origin of coordinates lies on dQ. Then 
for every f e L 2 (fl; l-R'l -1 ) there is a unique weak solution u G Hq (fi; K) to 
the Dirichlet problem A2.3\) , \2.4\ ) where K = x — y 2 . 

Proof. The proof follows the outline of the arguments in [17], Sec. 3. 
Defining a linear functional J/ by the formula 

J f (LO = (/, o , £ G C °° (O) , 

we estimate 

\Jf {Li) I < ||/|U2(n ; |K|-i)||ClU 2 (n ; |i<r|) < C||/||La(n;|K|-i)||-^||ir-i(n;JO> 

using Lemma 12.11 Thus Jf is a bounded linear functional on the subspace 
of if -1 (O; K) consisting of elements having the form L£ with £ e (fi) . 
Extending Jj to the closure of this subspace by Hahn-Banach arguments, 
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the Riesz Representation Theorem guarantees the existence of an element 
u G Hi (fl; K) for which 

(u,L£) = (f,0, 

where £ G Hq (fl; K) . There exists a unique, continuous, self-adjoint exten- 
sion L : Hq (fl; K) — > H -1 (fl; K) . As a result, if a sequence u n of smooth, 
compactly supported functions approximates u in the norm Hq (fl; K) , then 
Lu n converges in norm to an element / of H~ l (fl; K) . Taking the limit 

lim (u - u n , L£) — (f - /, f ) , 

we conclude that, because the left-hand side vanishes for all £ G Hi (fl; K) , 
the right-hand side must vanish as well. This proves the existence of a weak 
solution. Taking the difference of two weak solutions, we find that this differ- 
ence is zero in Hi (fl; K) by Lemma \2.1\ the linearity of L, and the weighted 
Poincare inequality [17J. This completes the proof of Theorem 12.21 

The unique-existence proofs of this section use estimates similar to those 
used in the proofs in Sec. 4 of [21] for existence alone. But the likelihood that 
the estimates would turn out to be similar appeared to be small on the basis 
of previous literature, and is rather surprising. In Sec. 5.1 of |24j it is shown 
that the estimates used to prove weak existence do not extend in an obvious 
way to proofs of uniqueness for the case K = x, in which the resonance 
curve is collinear with the flux line. Based on the physical discussion on 
p. 42 of (39], the collinear case would appear to be simpler than the case 
treated here, in which the two curves are tangent at an isolated point. But 
the simplicity of the case K = x arises largely from the fact that in this case 
the plasma behaves like a perpendicularly stratified medium, in which wave 
motion satisfies an ordinary differential equation; c.f. [T3]. The cautionary 
example of [21], Sec. 5.1, which suggests the difficulty of modifying the weak- 
existence methods to prove uniqueness in the case K = x, happens to fail 
in the case K = x — y 2 . As we have shown in this section, a modification 
of the weak existence estimates will in fact lead to a uniqueness proof for 
weak solutions to (12. 3p . ( 12. 4p for our choice of K. The case K = x remains 
interesting from several points of view, and we will return to its study in Sec. 
3.1. 

The restriction that the points of fl may not lie in the negative half-plane 
corresponds a requirement that boundary conditions be placed on one side of 
the flux line. As the resonance frequency is naturally restricted to the same 
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side of that line, the requirement seems to be compatible with the physical 
model. 

3 Strong solutions to open boundary value 
problems 

In this section we seek conditions sufficient for the existence of strong solu- 
tions in the cold plasma model. The conditions that we find will turn out to 
be extremely restrictive, but they are satisfied in a well known special case. 
Consider a system of the form 

Lu = f (3.1) 

for an unknown vector 

u = Oi (x,y),u 2 (x,y)) , 

and a known vector 

f = (/i (x,y),f 2 (?,y)) , 
where (x,y) G Q C M 2 . The operator L satisfies 

(Lu) 1 = K (x, y) u\ x + u 2y + zeroth-order terms, (3.2) 

(Lu) 2 = u X y - u 2x . (3.3) 

As in the preceding section, K (x, y) is continuously differentiable, negative 
on Q~ , positive on Q + , and zero on a parabolic region (the resonance curve) 
separating the elliptic and hyperbolic regions. If (/i, f 2 ) = (f, 0) , the compo- 
nents of the vector u are continuously differentiable, and u\ = u x , u 2 = u y for 
some twice-differentiable function u(x,y), then the first-order system (13. ip - 
( 13.31) reduces to a second-order scalar equation such as (12.31) . Because the 
emphasis in this section is on the form of the boundary conditions, the pres- 
ence or absence of zeroth-order terms will not affect the arguments provided 
the resulting system is symmetric positive. 

We say that a vector u = (ui,u 2 ) is in L 2 if each of its components is 
square-integrable. Such an object is a strong solution of an operator equation 
of the form (13. ip . with given boundary conditions, if there exists a sequence 
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u u of continuously differentiable vectors, satisfying the boundary conditions, 
for which \i u converges to u in L 2 and Lu u converges to f in L 2 . 

Sufficient conditions for a vector to be a strong solution were formulated 
by Friedrichs |8]. An operator L associated to an equation of the form 

Lu = A x vl x + A 2 u y + Bu, (3.4) 

where A 1 , A 2 , and B are matrices, is said to be symmetric positive if the 
matrices A 1 and A 2 are symmetric and the matrix 

Q = B*-l(Al + Al) 

is positive-definite, where B* is the symmetrization of the matrix B : 

B* = - (B + B T ) . 

The differential equation associated to a symmetric positive operator is also 
said to be symmetric positive. 

Boundary conditions for a symmetric positive equation can be given in 
terms of a matrix 

P = ntA\ dn + n 2 Af m , (3.5) 

where {n\,n2) are the components of the outward-pointing normal vector 
on dVt. The boundary is assumed to be twice- continuously differentiable. 
Denote by V the vector space identified with the range of u in the sense that, 
considered as a mapping, we have u : Q U dfl — > V. Let Af(x, y), (x, y) G dft, 
be a linear subspace of V and let Af(x, y) depend smoothly on x and y. A 
boundary condition u G is admissible if M is a maximal subspace of V with 
respect to non-negativity of the quadratic form (u, f3u) on the boundary. 

A set of sufficient conditions for admissibility is the existence of a decom- 
position ([Sj, Sec. 5) 

(3 = (3 + + (3_, (3.6) 

for which: the direct sum of the null spaces for /3 + and /3_ spans the restriction 
of V to the boundary; the ranges 9\± of /3± have only the vector u = in 
common; and the matrix \i = (3 + — (3- satisfies 

T 

H* = > 0. (3.7) 
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These conditions imply that the boundary condition 



/3_u = on ^ (3.8) 
is admissible for eq. ( 13. ip and the boundary condition 

w T /3l = 0ondn (3.9) 
is admissible for the adjoint problem 

L*w = g in fl. 

The linearity of the operator L and the admissibility conditions on the ma- 
trices (3± imply that both problems possess unique, strong solutions. 

Boundary conditions are semi-admissible if they satisfy properties (13.71) 
and (13. 8p . If f is in L 2 (Q) and (13 .ip is a symmetric positive equation having 
semi-admissible boundary conditions, then (13.11) possesses a weak solution in 
the ordinary sense: a vector u G L 2 (Q) such that 

/ (L*w)-udQ = I w-fdQ 
Jn Jn 

for all vectors w having continuously differentiable components and satisfying 
(PUD ([8J, Theorem 4.1). 

Writing the higher-order terms of eqs. (13. 2p . (13. 3p in the form 

M'^-OUH ;;)(:).■ (xio) 

we will derive admissible boundary conditions for the system (I3.1l) - ()3.3p . 

Slightly generalizing the type-change function of Sec. 2, we choose K (x, y) = 
x — a(y), where a(y) > is a continuously differentiable function of its ar- 
gument satisfying (c.f. |23j) 

tr(0) = ct'(O) = 0, (3.11) 
<r'(y)>0Vv>0, (3.12) 

and 

v'(y) < 0Vy < o. (3.13) 
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Let the operator L in (13.1 j) be given by 



(Lu) x = [x - er(y)] u Xx + w 2?/ + + k 2 u 2 , (3-14) 

(Lu) 2 = ui y - u 2x , (3.15) 

where K\ and «2 are constants. We note that equations of this kind satisfy 
all the physical requirements that originally led to the selection of eq. (II. ip 
in [39] as a model for cold plasma waves. 

By the elliptic portion dQ + of the boundary we mean points (x, y) of the 
domain boundary on which the type-change function K (x, y) is positive and 
by the hyperbolic portion dQ~ , boundary points for which the type-change 
function is negative. The sonic portion of the boundary consists of boundary 
points on which the type-change function vanishes. 

In this section we prove a revision and extension of [21], Theorem 9: 

Theorem 3.1. LetQ be a bounded, connected domain ofM 2 having C 2 bound- 
ary dfl. Let dQi be a (possibly empty and not necessarily proper) subset of 
dQ + . Let inequality A2.6\) be satisfied on dQ + \dQ± . On dQ± let 

b ni + cn 2 <0 (3.16) 

and on dQ\dQ + , let 

- &ri! + cn 2 > 0. (3.17) 

Let b(x, y) and c(x, y) satisfy 

b 2 + c 2 K^0 (3.18) 
on Q, with neither b nor c vanishing on Q + , and the inequalities: 

26ki - b x K -b + c y K - ca'(y) > in SI; (3.19) 

(26/«i — b x K — b + c y K — ca'(y)) (2ck 2 + b x — c y ) 

— (bn 2 + cki — c x K -c-b y f >0 inQ; (3.20) 

K (bm - cn 2 f + (cKm + bn 2 f < on dVt\dVL + . (3.21) 
Let L be given by (fjO^, ( Iff. 15\) . Let the Dirichlet condition 

— u\n 2 + u 2 n\ = (3.22) 
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be satisfied on dQ + \dfli and let the Neumann condition 

Ku\ti\ + u 2 n 2 = (3.23) 
be satisfied on dVt\. Then eqs. / TO) . ( fgTO ), HIM possess a strong solution 



on Q for every f G L 2 (Q). 

Proof. Multiply both sides of eqs. (13. ip . (I3.14p . (I3.15P by the matrix 

E = ( \ "f ) . (3,4) 

Condition ( I3.18P implies that E is invertible on Q, and conditions (I3.19p . 
(I3.20p imply that the resulting system is symmetric positive. 
For all points (x, y) G dQ, decompose the matrix 

oi~ ~\ _ ( K {bni - cn 2 ) cKni + bn 2 
I cKn\ + bn 2 — (brii — cn 2 ) 

into a matrix sum having the form = (3 + + 

On d{l + \d£l^ , decompose /3 into the submatrices 

/ Kbn x bn 2 

+ I Kcn\ cn 2 

and 

„ / —Kcn 2 Kcri\ 

\ bn 2 —bn\ 

Then /3_u = under boundary condition (13.221) . We have 



/i* = (bni + cn 2 ) 



K 
1 



so condition (I2.6P implies that the Dirichlet condition (I3.22j) is semi-admissible 

on dn\dn+. 

On dflf, choose 

-Kcn 2 Kcni 



1 bn 2 —bni 

and 

Kbni bn 2 



0- 



Kent cn 2 
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Then /3_u = under the Neumann boundary condition (I3.23p . and 

n ( K \ 
H = - (brix + cn 2 ) I g 1 J 

is positive semi-definite under condition (I3.16p . 

On dQ\dQ + , choose (3+ = (3 and take (3- to be the zero matrix. Then 
jj, = ji* = (3 and 

/in = K [bni - cn 2 ) . 

Because /in is non-negative by (13.171) . /i* is positive semi-definite by inequal- 
ity (13.211) . and no conditions need be imposed outside the elliptic portion of 
the boundary. 

This yields semi-admissibility. We now prove admissibility. 

On dQ + \dQf the null space of /3_ is composed of vectors satisfying the 
Dirichlet condition (13.221) . which is imposed on that boundary arc. The null 
space of (3 + is composed of vectors satisfying the adjoint condition (I3.23[) . 
On dQi , this relation is reversed. In order to show that the direct sum of 
these null spaces spans the two-dimensional space V\qq+, it is sufficient to 
show that the set 

\ W 1 

n 2 /ni J ' \ -Kn 1 /n 2 
is linearly independent there. Setting 

1 N \ ( 1 \ ( 



n 1 n 2 / ni +C2 \ -Km/n 2 J \ 



we find that C\ = —c 2 and 



_;<£±M ]=0 . (3.25) 



ni7i 2 

Equation (13.251) can only be satisfied on the elliptic boundary if c 2 = 0, 
implying that C\ = 0. Thus the direct sum of the null spaces of /3± on dfl + 
is linearly independent and must span V over that portion of the boundary. 

On dQ\dQ + , the null space of /3_ contains every 2- vector and the null 
space of (3 + contains only the zero vector; so on that boundary arc, their 
direct sum spans V. 

On dfl + \dflf , the range £H + of (3 + is the subset of the range of (3 for 
which 

v 2 rii — vin 2 = (3.26) 
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for (v i, v 2 ) G V; the range 9l_ of (3_ is the subset of 9^ for which 



Kv x n x + ^2^2 = 



(3.27) 



for (t>i,t>2) € V. Analogous assertions hold on dflf , in which the ranges of 
y{ + and are interchanged. Because if n\ and ri2 are not simultaneously 
zero the system (13.261) . (I3.27P has only the trivial solution t> 2 = v 1 = on 
dtt + , we conclude that 9t+ n 9t_ = {0} on <9fi + . 

On = {0}, so 9*+ n = {0} trivially. 

The invertibility of E under condition (13.181) completes the proof of The- 
orem EUl 

Remarks, i) By taking dflf to be either the empty set or all of dfl + , 
Theorem 3.1 implies the existence of strong solutions for either the open 
Dirichlet problem or the open Neumann problem. Because only the open 
Dirichlet problem was considered in Theorem 9 of [24] , Theorem 13.11 of this 
report extends that result to the open cases of the Neumann and mixed 
Dirichlet-Neumann problems. 

ii) A misprint in eq. (45) of [24J has been corrected in eq. ( 13.14]) . In 
Theorem 13. 1[ condition (I3.18P has been added to the list of hypotheses in 
Theorem 9 of [24] , the redundant condition (57) removed, and an error in eq. 
(59) corrected by eq. (I3.2ip of this report. The assumption that the boundary 
is piecewise smooth, which was default hypothesis in [24], seems to be too 
weak in general for strong solutions; see, however, [T2], [S], and [3T]-[3~4]. 

in) Only conditions (I3.19P and (I3.20P have anything to do with the cold 
plasma model. Otherwise, Theorem 13.11 is about interpreting Friedrichs' 
theory in the context of a collection of boundary arcs which are starlike 
with respect to a corresponding collection of vector fields. For example, the 
argument leading to eq. (I3.25P suggests that the Tricomi problem is strongly 
ill-posed under the hypotheses of the theorem, whatever the type-change 
function K. This is because in the Tricomi problem, data are given on both 
the elliptic boundary and a characteristic curve; but on characteristic curves, 
K satisfies 



Substituting this equation into eq. (13. 25 p . we find that the equation is satisfied 
on characteristic curves without requiring the constants C\ and C2 to be zero. 



K 




(3.28) 
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However, the theorem is less restrictive if the operator in (13. ip is given 

by 

{Lu) 1 = [x - cr(y)\ u lx - u 2y + Kim + k 2 u 2 , 

(Lu) 2 = -ui y + u 2x , (3.29) 

where, again, Ki and k 2 are constants. This variant also arises in the cold 
plasma model (see [22] and [30]) and is analogous to the variant of the Tricomi 
equation, 

t/Uxx ^yy 0) 

studied in various contexts by Friedrichs [§], Katsanis [llj, Sorokina [33J, 
[3lj . and Didenko [7]. In that case, choose 

*-(!?)■ 

Obvious modifications of conditions (I3.19P and (I3.20p guarantee that the 
equation 

ELu = Ef 

will be symmetric positive. Condition ( 13.1 8p must be replaced by the invert- 
ibility condition 

b 2 - c 2 K ± 0, 

which is restrictive on the subdomain Q + rather than on Q~ as in (13.181) . Most 
importantly, the discussion leading to Table 1 of [11] now applies, with only 
minor changes, and one can obtain a long list of possible starlike boundaries 
which result in strong solutions to suitably formulated problems of Dirichlet 
or Neumann type. In particular, one can formulate a Tricomi problem which 
is strongly well-posed. 

iv) The hypotheses of Theorem 13 .11 have a rather formal appearance. We 
expect them to be harsh, as the known singularity at the origin should dra- 
matically restrict the kinds of smoothness results that we can prove. But 
many of the conditions have natural interpretations. For example, inequali- 
ties (12.61) . (13.161) . and (I3.17P are satisfied whenever boundary arcs are starlike 
with respect to an appropriate vector field, and (I3.2ip is always satisfied on 
the characteristic boundary: 

Proposition 3.2. Let T be a characteristic curve for eq. Ii3. with the 
higher-order terms of the operator L satisfying Ii3.10\) . Then the left-hand 
side of inequality Ii3.21\) is identically zero on T. 
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Proof. We have, using eq. (I3.28p . 

(cKnx + bn 2 ) 2 = c 2 K 2 n\ + 2Kcbn x n 2 + b 2 n\ 



= -c 2 K 2 ^ + 2Kcbn x n 2 - b 2 Kn\ = -K (c 2 n 2 2 - 2cbn x n 2 + b 2 nj) 

= —K (cn 2 — brii) 2 . 

Substituting the extreme right-hand side of this equation into the second 
term of (13. 21 p completes the proof. 

3.1 An explicit example 

A simple example which illustrates the hypotheses of Theorem 13.11 can be 
constructed for the special case cr(y) = 0. In that case the system (I3.ip . 
(I3.14p . (13.15P can be reduced, by taking u\ = u x , u 2 = u y , and f = 0, to the 
Cinquini-Cibrario equation |5J 

xu xx + u yy + lower-order terms = 0. (3.30) 

In the context of the cold plasma model, this case corresponds to a resonance 
curve which is collinear with a flux line. In such situations eq. (II. ip can be 
replaced by an ordinary differential equation, as was discussed at the end of 
Sec. 2, so this choice is of little direct interest for the cold plasma model. But 
the Cinquini-Cibrario equation is interesting in its own right, in connection 
with normal forms for second-order linear elliptic- hyperbolic equations [3], 
[1]. Moreover, polar forms of eq. ( I3.30P arise in models of atmospheric and 
space plasmas — compare eq. (9) of [ID] , eq. (16) of [37], and eq. (19a) of 
[38] . with Sec. 3 of [DJ. (Such plasmas are cold, but in a relative rather than 
absolute sense.) 

In addition to taking a(y) = 0, choose K\ = 1; k 2 = 0; b = x + M, where 
M is a positive constant which is assumed to be large in comparison with 
all other parameters of the problem — in particular, & > OVi G fi; c = ey, 
where e is a small positive constant; (rii,n 2 ) = (dy/ds,—dx/ds), where s 
is arc length on the boundary. Inequalities (I3.18I) - (I3.20I) are satisfied for M 
sufficiently large. 
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Let the hyperbolic region il~ be bounded by intersecting characteristic 
curves originating on the sonic line. Condition ( I3.2ip is satisfied on dQ~ by 
Proposition 13.21 Condition (I3.17P is satisfied for M sufficiently large. As a 
concrete example, let dfl~ = T~ U T + , where 

r± = {(x,y) G Sl-\y = ±2 (yf=i -2)}. 

These curves intersect at the point (—4, 0) . Their intersection is not C 2 , but 
it can be easily "smoothed out" (by the addition of a small noncharacteristic 
curve connecting the points (— 4 + <5 , ±Si) for < S , 5\ « 1) without vio- 
lating either of the governing inequalities. Let the elliptic boundary dQ + be a 
smooth convex curve, symmetric about the with endpoints at (0, ±4) 

on the sonic line. Let the disconnected subset dQf of dQ + take the form 
of two small "smoothing curves," on which the slope of the tangent line to 
dfl + changes sign in order to prevent a cusp at the two endpoints. Inequality 
( 12. 6 p is satisfied on <9£7 + \fi^ and, again assuming that M is sufficiently large, 
inequality ( 13 . 1 6[) is satisfied on the two smoothing curves comprising dQ^. 

The domain Q of this construction is identical to the domain illustrated 
in Fig. 2 on p. 277 of j6], except that the cusps at points R, M, and N of 
that figure are smoothed out in Q near the points (—4, 0) and (0, ±4) . 

Theorem 13 . 1 1 implies that strong solutions to a homogeneous mixed Dirichlet- 
Neumann problem for the first-order inhomogeneous form of the Cinquini- 
Cibrario equation exist on this natural class of domains for any square- 
integrable forcing function. 

3.2 Symmetric positive operators on domains having 
multiply starlike boundaries 

On the basis of the observations in Remarks Hi) and iv), we reformulate 
Theorem 13.11 for an arbitrary type-change function — that is, for a smooth 
function K (x, y) which is positive on all points in a subset fl + of Q, negative 
on all points of a subset fl~ of fl, and zero on a smooth curve ^ G Q separating 
Q + and Q~ , where 

n = n + u or u & 

The proof of Theorem 13.11 will also prove: 

Corollary 3.3. Let Q be a bounded, connected domain of M 2 having C 2 
boundary dQ, oriented in a counterclockwise direction. Let dflf be a (possi- 
bly empty and not necessarily proper) subset of dfl + . Suppose that EL is a 



22 



symmetric positive operator, where L satisfies A3.1CA) (with the possible addi- 
tion of lower-order terms) and E satisfies Ij3.24\ ) with condition A3.18\) . Let 
dQ + \dQf be starlike with respect to the vector field V + = — (b(x, y), c(x, y)) ; 
let dfl^ be starlike with respect to the vector field V± = (b(x,y),c(x,y)) ; let 
dQ\dQ + be starlike with respect to the vector field V" = (b(x,y), —c(x,y)) . 
Let the union dQ\dQ + of the parabolic and hyperbolic boundaries be sub- 
characteristic in the sense of A3.21\) . Then the mixed boundary value prob- 
lem given by eqs. A3. 1\) - AH73)) . with condition A3. 22\) satisfied on dQ + \dQi 
and condition A3. 23\) satisfied on dQf, possesses a strong solution for every 

fei 2 (0). 
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